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Abstract 

Let ut = Uxx — q{x)ui < X < 1, t > 0, u{0, t) = 0, u{l, t) = a{t), u{x, 0) = 0, 
where a(t) is a given function vanishing for t > T, a{t) ^ 0, a{t)dt < oo. 
Suppose one measures the flux Ux{0,t) := 6o(t) for all t > 0. Does this information 
determine q{x) uniquely? Do the measurements of the flux Ux{l,t) := b{t) give 
more information about q{x) than 60 (^) does? 

The above questions are answered in this paper. 

1 Introduction 

Consider the heat transfer problem described by the equation 

Ut = Uxx — q{x)u, < x < 1, t > 0, (1.1) 

u{x,0) = 0, (1.2) 

w(0,t)=0, u{l,t)=a{t), (1.3) 

where a(t) is the prescribed temperature, and q{x) is a real-valued integrable function. 
Assume that a{t) is a pulse-type function, that is. 



a{t) =OfoTt>T, ! a{t)dt < 00, a{t) ^ 0. 
Jo 



[lA) 



In particular, one can choose a(t) to be the delta-function a{t) = S(t). Suppose one 
measures the flux at one of the ends of the rod, either measuring 

Ux{l,t):=b{t), (1.5) 



or 

UxiO,t) ■.= bo{t). 



(1.6) 



The questions that are answered in this paper are: 

1) Does the knowledge of a{t) and bo{t) for alH > determine q{x) uniquely? 

2) Does the knowledge of a{t) and b{t) for alH > determine q{x) uniquely? 

3) How does one calculate q{x) given a{t) and b(t)7 
The answers we give are: 

1) The knowledge of a{t) and bo{t) for all t > does not determine q{x) uniquely, 
in general. It does, if q{x) is symmetric with respect to the point x = |, that is, if 
q{x + |) = — x), or if q{x) is known on the interval [|, x]. 

2) The knowledge of a{t) and b(t) for alH > determines q{x) uniquely. 

3) An algorithm for computing q{x) given a{t) and b{t) is given. 

The answer to question 2) was given in [^] and earlier, under the additional assumption 
q{x) > 0, in [|I|. The answer to question 1) is new, as far as the author knows. An 
algorithm for computing q{x) is similar to the one described in 0]. 

In section 2 the answers to questions 1) and 2) are given, and an answer to question 
3) is given in section 3. 

2 Answer to questions 1) and 2). 

Let us Laplace-transform (1.1)-(1.3), (1.5) and (1.6). If v := v{x,X) := u{x,t)e~'^^dt, 
then 



where A{X), B{X) and Bo{X) are the Laplace transforms of a{t), b{t) and bo{t), respec- 
tively. 

Proposition 2.1. The data {^(A), -B(A)} known for a set of X > 0, which has a finite 
limiting point, determines q{x) uniquely. 

Proposition 2.2. The data {A{X) , Bq{X)} , known for all A > 0, does not determine 
q{x) uniquely, in general. 

If q{x) is known on the interval [|,l] then q{x) on the interval [O, |] is uniquely 
determined by the above data. 

Also, if q{\ — x) = g(| + x), then q{x) is uniquely determined on the interval [0, 1] by 
the above data. 

Proof of Proposition 2.1. Let (fi{x,iy) solve equation (2.1) with A = —u and satisfy the 
condition 



v" — q{x)v — Xv = 0, 0<a;<l 
v{0,X)=0, v{l,X) = A{X), 
v,{l,X) = B{X), 
t;,.(0,A) = So(A), 



(2.1) 
(2.2) 
(2.3) 
(2.4) 



ip{0,i^) = 0, <^'(0,z/) = l. 



(2.5) 
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The solution (/^(x, u) is an entire function of u and of A; = z/^/^ = iX^^'^ (see [Q). 
Since cp and v satisfy the first condition (2.2), one has: 

v{x,X) = c{X)^{x,-X), (2.6) 

where c(A) does not depend on x. Thus 

c(A)vp(l,-A) = A(A), c(A)v9'(l,-A)=S(A). (2.7) 

Note that t'(a;, A) may be not defined for some A, namely for some A, namely for —A = Aj, 
where Xj are the eigenvalues of the problem 

e^Pj:=ij'; + q{x)ij, = X,ijj, ^^.(0)=^,(1) = 0. (2.8) 

Since A > 0, the condition A = — Aj can be satisfies only if Xj < 0. There are at most 
finitely many negative eigenvalues of the selfadjoint Dirichlet operator i = — ^ + q{x) 
in H := L^[0, 1]. For the problem (2.1)-(2.2) to be solvable, when A = —Xj it is necessary 
and sufficient that the appropriate orthogonaltiy conditions are satisfied. Namely one 
finds 

^ A(X)Ml) 

v{x,X) = -Yl A + A ^^-^^ 

j=i 1 

For this series to be defined at A = —Xj > it is necessary and sufficient that 
A{—Xj) = 0. Note that "0^(1) 7^ by the uniqueness of the solution to the Cauchy 
problem (see (2.8)). 

Since we have assumed a{t) = for t > T, the function A{X) is an entire function 
of A on the complex A-plane. Therefore v{X) is well-defined as a meromorphic function 
of the parameter A with values in H. Note that problem (1. !)-(!. 3) is always solvable, 
but if the operator i has negative eigenvalues, then the solution to (1. !)-(!. 3) may grow 
exponentially as t — > +oo. 

From (2.7) one concludes 

B{X) ip'il - A) 



A{X) ^(1,-A)' 
since c(A) ^ 0. The zeros of the function 



(2.10) 



^{l,u)=0 (2.11) 
are precisely the Dirichlet eigenvalues Xj of i, while the zeros of the function 

ip'{l,iy) = (2.12) 
are precisely the eigenvalues of the problem 

iwj = fijWj, Wj{0) = 0, w'Jl) = 0. (2.13) 
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It is well known (see e.g. ||^) that the knowledge of {Xj} and {/ij} for all j determines 
q{x) uniquely because two spectra of i with the same homogeneous boundary condition 
at X = and two different homogeneous boundary condition at a: = 1, determine q{x) 
uniquely. 

The zeros of = ^^f^ are the numbers fij and only these numbers, while its poles 
are the numbers Xj and only these numbers. 
Proposition 2.1 is proved. 

□ 

Remark 2.1. A different proof of Proposition 2.1, based on Property C for ODE, is 
given in j^. 

Proof of Proposition 2.2. From (2.6) and (2.4) it follows that 

ciXMl, -A) = A{X), c(A)v5'(0, -A) = i?o(A). (2.14) 

Thus 

BojX) _ ^'iO,-X) _ 1 _ 1 

Aix) ipii,-x) y.(i,-A) ^ii,uy ■ ■ ^ ■ ' 

The poles of the function (2.15) are the eigenvalues Xj, and this is the only information 
one can get from (2.15). 

The knowledge of one spectrum {Xj} of i determines, roughly speaking, "half of the 
potential": namely, if q{x) is known on the interval [|,^], then the data {Xj} known 
for all j determine q{x) on [O, |] uniquely (see 0, [H, p). By the same reason if 
q{x + |) = ?(| — x) then q{x) is uniquely determined on [0, 1] by the set {Aj} known for 
allj. 

Proposition 2.2 is proved. □ 

The information in the data a{t) and bo{t) is equivalent to the information in the ratio 
^"^'^^ This is especially clear if one takes a{t) = 6{t) because in this case A(X) = 1 and 



A(A) 

Bp (A) _ 
A(A) <fi{l,u)- 



^sW. — — 1 — gQ f^i^g^^ t]2e information in the ratio is given just by one function (f{l,iy) 



Remark 2.2. In ^ and a general uniqueness result is obtained which says that if 
q{x) is known on [6, 1], < 6 < 1, where b is an arbitrary fixed number, then the 
set {Xm{j)} determines q{x) on [0,b] uniquely provided that a > 2b. Here Xm(j) is an 
arbitrary subset of {Xj} such that m{j) = ^(1 + Sj),Yl'jLi < oo. So, if m{j) = j, 
then a = 1, Sj = 0,b < ^. For & = | one gets the uniqueness result used in the proof of 
Proposition 2.2 and obtained in j^. 



Remark 2.3. From our arguments it follows that extra data (1.6) yields, roughly speak- 
ing, half of the information that data (1.5) yields, and therefore does not allow one to 
recover q{x) uniquely. 
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3 An algorithm for computing q{x) 



If {a{t), b(t)} are our data, one takes the Laplace transform and gets One calculates 
the zeros and poles of this function and gets the numbers {A-,} and {/ij}. In the literature 
(see 0) there is an algorithm for calculating the spectral function p(A) of the operator i 
from the knowledge of {Xj} U {/ij}- If p(A) is found, then the Gelfand-Levitan algorithm 
allows one to calculate q{x) from p(A). This algorithm is described in 0, [|], p. 

In this section we describe an algorithm which is a version of the one described in 
[5], pp. 297-299 (see [4], p. 57), which is quite different from the Gelfand-Levitan one and 
may be numerically more stable. 

Recall that 

r-^ sin(A^x) 
ip{x,u) = ipo{x,u) + K{x,y)ipo{y,'^)dy, ipo{x,i^) := ; , k = ^, (3.1) 



g(x) = 2^^^. (3.2) 



where K{x, y) is the transformation kernel, and 

dx 



Since Lp{y, Xj) = 0, one gets: 

Kil,y)^oiy,Xj)dy = -ipo{l,Xj), j = l,2,... (3.3) 



Since the set {(fioj} := {'^o{y, Aj)}vj, forms a Riesz basis of H = L'^[0, 1], relations (3.3) 
allow one to find K{l,y). 

Recall that a basis {hj} of a Hilbert space H is called a Riesz basis if there is a linear 
bounded map A and A"^ is a linear bounded operator on H, such that hj = Afj, where 
{fj} is an orthonormal basis of H (see p. 148). 

Numerically one may look for K{1, y) of the form 

J 

K(l,y) = J]c,(^o(y,A,), (3.4) 

substitute (3.4) into (3.3) and get a linear system for Cj,l < j < J. Here J is an arbitrary 
large positive integer. The matrix of the linear system is the Gram matrix 

sin{kjx) 



i^0j,^0m) ■= / ^Poj{x)(pornix)dx, ^Ojix) := ,kj = a/Aj, 

which is not ill-conditioned since {foj} forms a Riesz basis. 

Differentiate (3.1) with respect to x and set p = fij x = 1 to get 

= v?[j(l,/x,) + K(l,l)^o(l,Aii)+ /" K,{l,y)^o,{y)dy, j = l,2,.... (3.5) 
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These equations determine uniquely Kx{l,y), since v?q(1, /ij), 1) and Lpo{l,fij) are 
known numbers. Thus we can compute K{l,y) and Kx{l,y), < y < 1, from the data 
{a{t),b{t)}. 

If K{l,t) and Kx{l,t) are known, then one can derive a Volterra integral equation 
for the unknown U := {q{x), K{x,y)} (see @], p. 56, and 0). 

In it is proved that this equation can be solved by iterations, and therefore q{x) 
can be computed by an iterative process. 

For convenience of the reader we write down the above integral equation for U := 
{q{x), K{x, y)} and an iterative process for the solution of this equation: 



U = W{U) + h, 



where 



W{U) :-- 



-2 f q{s)K{s,2x - s)ds 

J X 



\ 2 



q{s)K{s, t)dsdt 



(3.6) 



(3.7) 



J 



Dxy is the region bounded by the straight lines s = l^t — y = s — x, and t — y = x — s 
on the (s, t) plane, 

h = \9 \ , (31 



fix) := 2[Ky{l, 2x~l) + K,{1, 2x - 1)], 

( ^ Kil,y + x~l) + Kil,y-x + l) 1 fy~^+' 
9{x,y) = / Ks{l,t)dt. 

^ ^ Jy+x~l 



(3.9) 
(3.10) 



Note that / and g are computable from the data K{l,x) and Xa;(l,x), and (3.6) is a 
nonlinear Volterra-type equation for the unknown q{x) and K{x,y). 
It is proved in [Q, p 57 (and in 0) that the iterative process 



Un+l = W{Un) +h, Uo = h 

q{x) 



converges (at a rate of geometric series) to U (x) 



K{x,y) 

the functional space in the norm of which the convergence hold are given in 



(3.11) 

. The details concerning 



References 

[1] Denisov, A., Introduction to the theory of inverse problems, MGU Press, Moscow, 
1994. 

[2] Hochstadt, H., Lieberman, B., An inverse Sturm-Liouville problem with mixed data, 
SIAM J. of Appl. Math., 34, 1976, 676-680. 



6 



[3] Levitan, B.M., Inverse Strum- Liouville problems, VNU Press, Utrecht, The Nether- 
lands, 1987. 



[4] Ramm, A.G., Property C for ODE and applications to inverse problems, in the 
book "Operator Theory and Its Applications", Amer. Math. Soc, Fields Institute 
Communications vol. 25, (2000), pp. 15- 75, Providence, RI. (editors A. G. Ramm, 
P.N.Shivakumar, A.V.Strauss). 

[5] Ramm, A.G., Multidimensional inverse scattering problems, Long- 
man/Wiley, New York, 1992, pp. 1-385. 

[6] Ramm, A.G., Property C for ODE and applications to inverse scattering, Zeit. fiir 
Angew. Analysis, 18, N2, (1999), 331-348. 

[7] Ramm, A.G., Recovery of compactly supported spherically symmetric potentials 
from the phase shift of s-wave. In the book: Spectral and scattering theory. Plenum 
publishers. New York, 1998 (ed. A.G.Ramm), pp.111-130. 

[8] Ramm, A.G., Scattering by obstacles. Reidel, Dordrecht, 1986, pp. 1-442. 



7 



